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Crystallographic orientations determine the optical, electrical, mechanical, and thermal properties of 
crystals. Control of crystallographic orientations has been studied by changing the growth parameters, 
including temperature, pressure, proportion of precursors, and surface conditions. However, molecular 
dynamic mechanisms underlying these controls remain largely unknown. Here we achieved control of 
crystallographic orientations in diamond growth through a joint experimental and theoretical study of laser 
resonant vibrational excitation of precursor molecules (ethylene). Resonant vibrational excitation of the 
ethylene molecules using a wavelength-tunable C0 2 laser steers the chemical reactions and promotes 
proportion of intermediate oxide species, which results in preferential growth of {100}-oriented diamond 
films and diamond single crystals in open air. Quantum molecular dynamic simulations and calculations of 
chemisorption energies of radicals detected from our mass-spectroscopy experiment provide an in-depth 
understanding of molecular reaction mechanisms in the steering of chemical reactions and control of 
crystallographic orientations. This finding opens up a new avenue for controlled chemical vapor deposition 
of crystals through resonant vibrational excitations to steer surface chemistry. 



Laser control of chemical reactions through resonant excitations has been extensively studied in bimolecular 
collisions, unimolecular decompositions, and gas-solid reactions 15 . However, the effects of laser control of 
reactions are usually too weak to be significant for material synthesis 6 . Control of chemical reactions via 
resonant vibrational excitation of precursor molecules is promising in that molecular vibrations are directly 
related to the rearrangement of their bonds to become products 1,5 . Lasers are ideal energy sources for vibrational 
excitations because their narrow bandwidth prepares a single vibrational eigenstate 5 , which makes precise control 
of reaction pathways feasible. One approach to laser control is active intervention, in which the phases of reagent 
motions are controlled during the course of the reaction 6 . This active intervention may preferentially lead 
reagents to one of many different reaction routes. In this study, we introduced active intervention to material 
synthesis, and achieved control of chemical reaction and hence controlled crystallographic orientations through 
laser resonant vibrational excitations of precursor molecules in chemical vapor deposition (CVD) of diamond. 

Due to its extreme properties, diamond has found widespread applications, including heat spreaders, optical 
windows, biosensors, and nano- or microelectromechanical systems (N/MEMS) 71 ", as well as its promising 
potentials in quantum information technology 1113 . Many properties of diamond, including thermal conductivity, 
optical behavior, and mechanical properties, depend on diamond crystallographic orientations 1415 . Preferential 
growth of { 100}-textured diamond films has been pursued for years because of their superior properties over those 
of the {111}- and {110}-textured diamonds 1623 . Mechanically, {100}-textured diamond fdms have lower rough- 
ness and higher wear resistance as compared with other crystallographic directions 16 . Optical properties of {100}- 
textured films are superior to those of { 1 1 1 }-textured films in terms of refractive index and extinction coefficient 16 . 
Research efforts have been exerted to achieve orientation-controlled growth of diamond. Liu et al. have studied 
control of diamond textures by varying the growth parameters, including substrate temperature, precursor 
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Figure 1 | Absorption of laser energy by the C2H4/C2H2/O2 flame, (a), Schematic of the experimental setup, (b), Absorption of laser energy by the 
C2H4/C2H2/O2 flame as a function of laser wavelength in the C0 2 laser 00°1-10°0 emission band. The error bars indicate standard deviations, 
(c), Vibration-rotation transitions of an ethylene molecule with C0 2 laser excitations at different wavelengths, indicating that the C0 2 laser at both 10.22 
and 10.532 urn excites the molecules to a higher vibrational level, whereas the 10.22 um excitation provides additional rotational energy in the molecules. 



concentration, and films thickness 14 . By applying bias-enhanced 
nucleation 20 ' 21 and introducing nitrogen into the growth precur- 
sors 2226 , {100}-textured diamond films have also been successfully 
deposited. These methods are capable of producing {100} -textured 
diamond films with good reproducibility. However, the molecular 
reaction mechanism remains elusive and further investigations are 
needed to obtain in-depth understanding on the control of crystal- 
lographic orientation. In our previous study, resonant excitations of 
precursor molecules have been applied in the synthesis of diamond 
films and crystals 27,28 . Diamond growth rate and quality were 
obviously promoted. Diamond crystals with sizes up to 5 mm were 
successfully obtained in open air 28 . In-depth understanding of the 
laser control of chemical reactions and control of diamond crystal- 
lographic orientations are investigated in this research. 

Results and discussion 

We introduced resonant excitations of precursor molecules using a 
wavelength-tunable C0 2 laser to steer the chemical reactions in dia- 
mond formation which in turn resulted in preferential growth of 
{100} -textured diamond of high phase purity. The diamond growth 
process was steered by resonantly exciting ethylene (C 2 H 4 ) molecules 
to higher vibrational and rotational states. A schematic illustration of 
the experimental setup is shown in Figure la. The detailed experi- 
mental method is described in the Methods section. The wavelength 
of the tunable C0 2 laser ranges from 9.2 to 10.9 (.im. Figure lb shows 



a representative absorption spectrum of the laser energy by the C 2 H 4 / 
C 2 H 2 /0 2 flame with respect to the laser wavelength in the C0 2 laser 
00°1-10°0 emission band. Several absorption peaks were observed at 
different wavelengths. The strongest peak at 10.532 urn corresponds 
to the Q branch (AJ = 0) of a fundamental vibration mode (v 7 , CH 2 - 
wagging) of the ethylene molecules 25 . Another absorption peak at 
10.22 (im corresponds to the R branch (AJ = 1) of the CH 2 -wagging 
mode 29,30 , which means that the molecules are excited to a higher 
level of vibrational state with a higher level of rotation energy, as 
schematically indicated in Figure lc. 

We compared the excitations of ethylene molecules in the com- 
bustion flame by laser irradiation at 10.22, 10.333 (non-resonant 
wavelength), and 10.532 |j.m, and deposited diamond films with all 
other parameters the same. Figure 2a-d show representative scan- 
ning electron microscopic (SEM) images of diamond films deposited 
for 15 (upper) and 60 min. (lower), respectively. Diamond films 
deposited without laser excitation were also included for compar- 
ison. Thicknesses of these diamond films are given in Figure S 1 of the 
Supporting Information. It was observed that diamond films depos- 
ited by the 10.22 u.m laser excitation exhibit predominantly {100} 
surfaces, while those deposited under the other three conditions 
showed random orientations. This phenomenon does not change 
with deposition time, as demonstrated by similar film morphology 
under each condition of different deposition time (15 and 60 min.). 
It is demonstrated that excitation of ethylene molecules through the 
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Figure 2 | SEM images of the diamond morphologies, (a-d), Diamond films deposited without laser (a), with C0 2 laser irradiation at 10.22 (b), 10.333 
(c), and 10.532 urn (d). (upper: 15-min. deposition, lower: 60-min. deposition), (e), (f), Single-crystal diamond grown for 5 hr. with C0 2 laser 
excitation at 10.22 um. 



vibrational excitation of ethylene molecules at 10.22 (tm leads to a 
preferential growth of {100} -textured diamond films. 

Under the resonant excitation at 10.22 um, single crystal diamond 
with a size of 300 (tm was grown in 5 hr. in open air, as shown in 
Figure 2e. Different from the diamond crystals deposited through 
vibrational excitation of ethylene at 10.532 um 28 , this diamond crys- 
tal exhibits high uniformity. A flat {100} surface as shown in Figure 2f 
demonstrated the preferential growth of the {100} surface for long 
time deposition. No twinning was observable in the SEM images on 
the {100} surface. Twinnings abundantly exist on the {111} surfaces 
(the four surfaces surrounding the {100} surface), which is believed to 
be caused by the faster growth in the (111) compared to that in the 
(100) direction 15 . The crystal shape of diamond facets is determined 
by the growth rate of the (100) direction (V 10 o) and that of the (111) 
direction (Vm), as indicated by the diamond growth parameter, 

« = \/3 — — 14,32 . The facets that show up at the final stage of crystal 
Vm 

growth are those along which the growth velocity is the least. The 
growth velocity of the (111) direction should be greater than that of 
the (100) direction in order to obtain the predominantly {100} sur- 
faces. It has been widely accepted in diamond growth that higher 
growth velocity in certain sector will result in lower phase purity in 
that sector, which results in the twinnings on the {111} surfaces 15,31,32 . 
A lower concentration of twinnings on the {100} surface indicates 
high quality of the {100} -oriented diamond crystals. 

Raman spectroscopy was performed to determine the phase purity 
of the deposited diamond films and crystals. Figure 3a shows Raman 
spectra of the diamond films deposited for 60 min. under different 
excitation conditions. The shift of the diamond peak from 1332 
(phase-pure natural diamond) to 1335 cm 1 is ascribed to the resid- 
ual stresses in the diamond films caused by differences in the coeffi- 
cients of thermal expansion (CTE) between the diamond films and 
the substrates. The D-band is related to disordered carbon structures, 
and G-band corresponds to sp 2 graphite phases in the diamond films. 
It is obvious that both D- and G-bands are suppressed in diamond 
films deposited with laser excitations at 10.22 and 10.532 (tm, indi- 
cating an increased purity of diamond sp 3 bonding. The spectrum of 
the diamond film deposited with laser excitations at 10.22 um 



exhibits even lower D- and G- bands, indicating that the {100} sur- 
faces have better phase purity than other surfaces. This is in agree- 
ment with the transmission electron microscopic (TEM) study of 
different diamond growth sectors, which revealed that the (100) 
sector has low density of dislocations, while the (111) sectors are 
highly defective, containing a high density of dislocations and stack- 
ing faults 15 . Figure 3b shows the Raman spectra of diamond crystals, 
including those deposited with 10.22 (shown in Figure 2e) and 
10.532 (tm excitations (reported in Ref. 25) as well as a type la natural 
diamond. The diamond peaks of these three crystals are all located at 
1332 cm" 1 . The D- and G-bands are invisible in all spectra, which 
indicate very high purity of the diamond sp 3 bonding. The full width 
at half maximum (FWHM) value of the {100}-oriented diamond is 
3.3 cm -1 , which is similar to that of the natural diamond (3.1 cm -1 ), 
much narrower than that of the diamond crystal deposited with the 
10.532 (tm excitation (4.5 cm -1 ). Therefore, the {100}-oriented dia- 
mond crystal deposited with vibrational excitation of ethylene at 
10.22 um has high phase purity and high crystal quality. 

Mass spectroscopy (MS) of the C2H4/C2H2/O2 combustion flame 
used for diamond growth was performed under different excitation 
conditions to investigate mechanisms of the preferential growth of 
{100}-textured diamond, with more details described in the Methods 
section. Ionization of species occurs in the combustion flame, mak- 
ing the flame suitable for direct analysis using MS. Figure 4a shows a 
representative mass spectrum of the C2H4/C2H2/O2 flame, in which 
positive ions with m/z values ranging from 12 to 47 are detected and 
assigned as CH X + (x = 0 ~ 3), C 2 H X + (x = 0 ~ 5), C 3 H X + (x = 0 ~ 5), 
CH x O + (x =1 — 3) and C 2 H x O + (x = 2 — 5) ions. Considering 
differences in ionization potentials among different species, the 
height of each line does not necessarily represent the absolute con- 
centration of each species. Nevertheless, MS analysis of the flame is 
valuable to determine the evolvement of species concentration under 
different excitation conditions, given that the variation of ions is 
proportional to that of related species. Figure 4b shows a chromato- 
gram of the total ion current (TIC) of the C2H4/C2H2/O2 flame under 
different excitation conditions, which provides the intensity of all 
detected ions in the flame. The TIC of the flame is a summation of 
all detected ions in the flame. By using the analysis tools as described 
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Figure 3 | Raman spectra of the diamond samples, (a), Diamond films deposited for 60 min without laser, with C0 2 laser excitation at 10.22, 10.333 
and 10.532 urn. (b), Diamond crystals grown for 5 hr with C0 2 laser excitation at 10.22 and 10.532 |im. Natural diamond (type la) is also 
characterized as comparison. Insets respectively show the FWHM of the diamond peak for each crystal. 



in the Methods section, the chromatogram of each detected ion can 
be derived from the TIC signal. The chromatogram of each ion has 
similar shape with that of the TIC. Figure S3, given in the 
Supplementary Information, provides chromatograms of CH 3 + , 
CH 3 0 + , C 2 H 3 + , and C 2 H 3 0 + as examples. The chromatograms of 
other species (CH, CH 2 , C 2 H, C 2 H 2 , CH 2 0, C 2 H 2 0, etc.) have sim- 
ilar shapes. The intensity of each ion can be obtained by calculating 
the average value of the corresponding scans. The concentration of 
each ion was derived by dividing its intensity with the TIC value in 
order to eliminate the variation of its absolute value in different 
measurements. Because species containing one or two carbons are 
of interest to diamond growth 32 , relative concentrations of CH X + 
(x = 0 ~ 3), C 2 H X + (x = 0 ~ 5), CH x O + (x = 1 ~ 3), and 
C 2 H x O + (x = 2 ~ 5) ions were calculated and summarized respect- 
ively, as shown in Figures 4c and 4d. The relatively large deviation of 
the species concentration is mainly caused by flame fluctuation dur- 
ing the mass spectrometry (MS) measurement. The uncertainty 
caused by the deviation is suppressed in two aspects. Firstly, mass 
flow controller (MFC) with more precise control was used, which 
suppressed the flame fluctuation. Secondly, number of the measure- 
ments (more than twenty) was increased for each wavelength. The 
mean value for each wavelength is hence believed to be representative 
for the trends, although relatively large deviation exists. 

In Figure 4c, the concentrations of both CH X and C 2 H X species 
were increased by laser excitations at 10.22 and 10.532 um. The 
CH x O, and C 2 H x O species show similar trends, as shown in 
Figure 4d. However, a detailed comparison between Figure 4c and 
4d reveals a slight difference between the species intensities with laser 
excitations at 10.22 and 10.532 um. Concentrations of the CH X and 
C 2 H X species are the highest at 10.532 um, whereas those of the 
CH x O and C 2 H x O species are the highest at 10.22 um. As discussed 
in the beginning of this section, the excitation of ethylene molecules 
by the C0 2 laser at 10.22 fim corresponds to the R branch (AJ = l)of 
the CH 2 -wagging mode, which means that the molecules are excited 
to a higher level of vibrational state with a higher level of rotation 



energy. Because the energy transfer by the rotational-translational 
channel through binary molecular collision is much faster than that 
by the rotational-vibrational channel 33 , higher rotational energy will 
result in higher translational energy, which is reflected by higher 
temperature. That means the excitation of ethylene molecules by 
the C0 2 laser at 10.22 um resulted in higher flame temperature than 
the excitations at other wavelengths (with same quantity of laser 
energy absorbed). The increased flame temperature consequently 
resulted in promoted oxidation of the carbon hydride species. In 
other words, concentration of the oxygen- containing species 
(CH x O, C 2 H x O, etc.) was increased by higher flame temperature 
caused by the 10.22 um C0 2 laser excitation. This explanation has 
been confirmed by quantum molecular dynamics (QMD) simula- 
tions for exploration of intermediate species in the C 2 H 4 /C 2 H 2 /0 2 
flames, shown in Movie SI and summarized in Figure S4. It is indi- 
cated in the simulation that as the temperature of the flame is 
increased, the CH x O (x = 1 ~ 3) and C 2 H x O (x = 2 ~ 5) species 
appear with larger population through the combustion process. 
Details of the QMD simulation are given in the Supplementary 
Information. 

To gain molecular insights into relative abundances of produced 
species on different diamond surface, we computed the chemisorp- 
tion energies of all species detected in the mass spectrometry experi- 
ment (CH X (x = 0 ~ 3), C 2 H X (x = 0 ~ 5), CH x O (x = 1 ~ 3), and 
C 2 H x O (x = 2 ~ 5)) on the diamond C{100} and C{111} surface, 
respectively. Periodic density-functional theory (DFT) methods are 
employed (see Methods section). The chemisorption energies are 
shown in Table 1. Note that the chemisorption capacity can be used 
as an indicator of surface reactivity 34,35 . Stronger chemisorptions lead 
to higher surface reaction activity. As shown from Table 1, the 
chemisorption of all species (except C 2 H 5 0 with the least population, 
which is due to the perfect bond-length match with the C{100}-2 X 1 
reconstruction surface) on the C{111} surface is much stronger 
than on the C{100} surface, indicating higher reactivity on C{111} 
surface and hence higher growth in the (111) direction than that in 
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Figure 4 | Mass spectrometry of the C2H4/C2H2/O2 flame, (a), Mass spectrum of the flame irradiated with C0 2 laser at 10.22 ^m. (b), Total ion 
current of the flame under different excitation conditions, (c), (d), Relative concentrations of CH X (x = 0 ~ 3), C 2 H X (x = 0 ~ 5) ions (c), and CH x O 
(x = 1 ~ 3), C 2 H x O (x = 2 ~ 5) ions (d). 



the (100) direction. This favorable adsorption of hydrocarbon species 
on C{111}-1 X 1 surface over the C{100}-2 X 1 reconstruction sur- 
face is very similar to H- and OH- adsorptions on the relevant dia- 
mond surfaces 36 3S . Keep in mind that faster growth of one axial 
direction of crystals will result in larger areas of adjacent facets 14,32 . 
In this case, the growth rate in the (1 1 1) direction is higher than that 
in the (100) direction, resulting in larger area of the adjacent C{100} 
facets, which explains the predominant appearance of the diamond 
{100} surfaces. 

To summarize, we have demonstrated a laser-assisted combustion 
process to synthesize {100} -oriented diamond films and single crys- 
tals through resonant vibrational excitation of precursor molecules. 
The C0 2 laser beam at 10.22 um excites the ethylene molecules to a 
higher vibrational state with higher rotational energy in the combus- 
tion flame, which increases the relative concentrations of the CH x O 
and C 2 H x O species. Theoretical simulations of the reaction pathway 
indicate that these species react more easily with diamond {111} 
surfaces than with {100} surfaces, leading to a fast growth in the 
(111) direction and resulting in {100}-oriented diamond. This find- 
ing opens up a new avenue for controlled chemical vapor deposition 
of crystals through resonant vibrational excitation of precursor 
molecules. 



Table 1 | Calculated chemisorption energies (eV) on Diamond 
{1 00} and {1 1 1 } surface with PBE/DNP method 



. . Chemisorption energy (eV) 
lecule 





On C{100} 


On C{111} 


CHO 


-2.58 


-3.28 


CH z O 


-2.64 


-3.27 


CH3O 


-2.83 


-3.53 


C 2 H 


-4.43 


-5.53 


C 2 H 2 


-3.75 


-4.60 


C 2 H3 


-3.84 


-5.78 


C 2 H 2 0 


-3.13 


-3.99 


C 2 H 3 0 


-4.69 


-6.49 


C 2 H 5 0 


-6.09 


-3.53 


C 3 H 


-4.04 


-4.34 


C3H 2 


-4.06 


-4.93 


C3H3 


-2.35 


-2.76 


C3H4 


-3.62 


-4.34 


C3H5 


-5.18 


-6.35 



SCIENTIFIC REPORTS | 4:4581 | DOI: 10.1038/srep04581 



5 



Methods 

Diamond synthesis. A combustion torch with a 1.5 mm orifice tip was used to 
produce the flame. The precursor gases consisted of a mixture of C 2 H 4 (99.999%), 
C 2 H 2 (99.6%), and 0 2 (99.996%), which were mixed in the torch through three gas 
flow meters (B7920V, Spec- Air Gases & Technologies). The gas flow rates of C 2 H 4 , 
C 2 H 2 , and 0 2 were 0.62, 0.62, and 1 .20 standard liters per minute (slpm), respectively. 
The C0 2 laser beam used to irradiate the flame was parallel with the substrate 
surfaces. A ZnSe convex lens with a focal length of 256 mm was used to focus the laser 
beam from its original diameter of — 13 mm to ~-2 mm which is similar to the 
diameter of the flame {Figure la). The laser power density at the flame was around 
3 X 10 4 W/cm 2 . A tungsten carbide (WC) substrate (BS-6S, Basic Carbide Corp., 
containing 6% cobalt) with a dimension of 12.7 X 12.7 X 1.6 mm 3 was placed on a 
water-cooled brass box mounted on an X-Y-Z stage. The surface roughness of the WC 
substrate was 400 nm. The substrate was cleaned in a supersonic bath of acetone for 
15 min before diamond deposition. The distance between the substrate surface and 
the inner cone of the flames was about 0.5 mm in all experiments. The temperature of 
the substrate during deposition was maintained at 770 ~ 780 J C and monitored by a 
noncontact pyrometer (OS3752, Omega Engineering, Inc.). Raman spectroscopy was 
performed for the diamond samples using a Renishaw Raman spectrometer (in Via 
9F5803) with an argon ion laser as the excitation source (514.5 nm, 100 mW). 

Mass spectrometry. Positive ions in the C 2 H 4 /C 2 H 2 /0 2 flame were detected using a 
time- of- flight mass spectrometer (MS, AccuTOF™, JEOL USA, Inc.), as shown in 
Figure S2 in the Supplementary Information. A stainless steel orifice with an inner 
diameter of 400 urn on the mass spectrometer collects ions in open air. The 
combustion torch was fixed on a motorized X-Y-Z stage and the relative position of 
the flame to the orifice was precisely adjusted so that the tip of the flame inner cone 
was right at front of the orifice under all excitation conditions. The flame-to- orifice 
distance for the MS detection was the same as the flame- to -substrate distance for 
diamond deposition. In this way, the detected species are believed to be the same as 
what exist in the diamond growth. The MS data were analyzed using the TSSPro 
software (Shrader Analytical and Consulting Laboratories, Inc. Version 3.0) and the 
MS Tools software (JEOL USA, Inc.). 

Computation of chemisorption. The chemisorption energies were calculated using 
the density functional theory (DFT). We constructed a slab containing 8 atomic layers 
of 4 X 4 diamond {111}-1 X 1 surface (where the lower 4 atomic layers are fixed and 
the top 4 atomic layers are fully relaxed) and another slab containing 9 atomic layers 
of 4 X 4 {100}-2 X 1 reconstruction surface (where the lower 4 atomic layers are fixed 
and top 5 atomic layers are fully relaxed). Each layer has 16 carbon atoms and the 
lower 4 layers are not terminated. Note that the slab models constructed here are very 
close to the models reported in the literatures 36-38 . The vacuum slab in the supercell 
was set as 20 A. The DFT calculations are based on the Perdew-Burke-Ernzerhof 
(PBE) 39 functional at double numerical polarized (DNP) basis sets 40 41 . The DNP basis 
sets are considered as superior as the Gaussian 6-311 + G(3 df,2 pd) basis sets 42 , 
confirming reliability of our DFT calculations. The energy convergence was set as 
10~ 6 Hartree, and the energy convergence for the geometrical optimization was set as 
1 X 10~ 5 Hartree. The chemisorption energies were evaluated as: E(chemisorption) = 
E(C n H m Ok — surface) — E(C n H m Ok) — E(surface). All calculations were carried out 
using the DMol 3 6.0 program software 40 ' 41 . 
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